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ABSTRACT 

In the present paper we develop an algorithm allowing to calculate line-of-sight ve- 
locity dispersions in an axisymmetric galaxy outside of the galactic plane. When con- 
structing a self-consistent model, we take into account the galactic surface brightness 
distribution, stellar rotation curve and velocity dispersions. We assume that the ve- 
locity dispersion ellipsoid is triaxial and lies under a certain angle with respect to the 
galactic plane. This algorithm is applied to a Sa galaxy NGC 4594 — M 104, for which 
there exist velocity dispersion measurements outside of the galactic major axis. The 
mass distribution model is constructed in two stages. In the first stage we construct 
a luminosity distribution model, where only galactic surface brightness distribution is 
taken into account. Here we assume the galaxy to consist of the nucleus, the bulge, the 
disc and the stellar metal-poor halo and determine structure parameters of these com- 
ponents. Thereafter, in the second stage we develop on the basis of the Jeans equations 
a detailed mass distribution model and calculate line-of-sight velocity dispersions and 
the stellar rotation curve. Here a dark matter halo is added to visible components. 
Calculated dispersions are compared with observations along different slit positions 
perpendicular and parallel to the galactic major axis. In the best-fitting model velocity 
dispersion ellipsoids are radially elongated with ae/an ~ 0.9 — 0.4, (Jz/ctr — 0.7 — 0.4, 
and lie under the angles ^ 30° with respect to the galactic equatorial plane. Outside 
the galactic plane velocity dispersion behaviour is more sensitive to the dark matter 
density distribution and allows to estimate dark halo parameters. For visible matter 
the total M/Lb = 4.5 ±1.2, M/Lr = 3.1 ±0.7. The central density of the dark matter 
halo is /odm(O) = 0.033 Mqpc^^ 

Key words: galaxies: individual: NGC 4594 - galaxies: kinematics and dynamics - 
galaxies: spiral - galaxies: structure - cosmology: dark matter. 



1 INTRODUCTION 

The study of the dark matter (DM) halo den- 
sity distribution aUows us to constrain possible 
galaxy formation models and large scale struc- 
ture formation scenarios jNa varro fc Stcinmctz 2000|; 
i Khairul Alam. BuUock fc Weinbera .2002,: Gentile et all 
12004^ . For this kind of analysis, it is necessary to know 
both the distribution of visible and dark matter. Without 
additional assumptions rotation curve data alone are not 
sufficient to di s crimin ate between these two kinds of matter 
jPutton et alJ |200^. It does not suffice either to use 
additionally velocity dispersions along the major axis. 

Realistic mass and light distribution models must be 
consistent, i.e. the same model must describe the luminos- 
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ity distribution and kinematics. Three main classes of self- 
consistent mass distribution models can be discriminated: 
the Jeans equations based models, the specific phase space 
density distribution models and the Schwarzschild orbit su- 
perposition based models. 

Mass distribution models based on solving the Jeans 
equations have an advantage that the equations contain ex- 
plicitly observed functions - velocity dispersions. On the 
other hand, there are three equations, but at least five un- 
known functions (three dispersion components, centroid ve- 
locity and the velocity dispersion ellipsoid orientation pa- 
rameter) and thus the system of equations is not closed. 
In addition, the use of the Jeans equations neglects possi- 
ble deviations of velocity distributions from Gaussians and 
does not garantee that the derived dynamical model has 
non-negative phase density distribution everywhere. How- 
ever, within certain approximations the Jeans equations 
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are widely used for the construction of mass distribution 
models. In the case of spherical systems with biaxial ve- 
locity dispersion ellipsoi ds, such models have b e en con- 
struct e d, for exa mple, by'Binne^^^M^mo^jEls^),!^^^!^ 
|l983),|Gerha3 Cl991), Tr emaine et al.1 Jl994 ."ln the case 
of flattened systems with biaxial velocity dispersion ellip- 
soids, a general algorithm f or the solution of the Jeans 
equations was developed bv [ B inneyi Davies fc lUinerworthl 
l)l99(tl . lcinzano fc van der MareJ (|l994 ). Another algorithm 
in the context of the multi- Gaussian expan si on (M GE) 
formalism was developed by lEmsellem et al.l il994l) . An 
approximation for cool stellar discs (random motions are 
small when compared with ro tation) has been developed by 
lAmendt fc CudderfordI Jl99ll) . 

Dynamical models with a specific phase density 
distribution have the advantage that velocity dispersion 
anisotropy can be calculated directly. On the other hand, 
due to rather complicated analytical calculations, only 
rather limited classes of distribution functions can be stud- 
ied . Spherical models of this kind ha v e bee n constructed 
bvlCaroUo. de Zeeuw fc van der Marell lll995ll . iBertin et alJ 
In the case of an axisymmetric density distribution, 
velocity dispersion profiles have been calculated for certain 
specific mass and phase densi ty distr ib ution f or ms b y 
van der Marel. Binnev fc DaviesI Jl99d). lEvansI Jl993l) 
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1996[) . lde Zeeuw. Evans fc Schwarzschildl Jl99d) . iMerritt 
19961) , I An fc Evand fcOOd) and others. A special 



case IS an 



analytical solution with three integrals of motion for some 
specific potentials: an axisymmetric model with a potential 
in the Sta ckel form jPej onghe fc de Zeeuwlll988f) . isochrone 
potential JPehnen fc Gerhardll993l) . 

Probably the most complete class of dynamical mod- 
els has been developed on t he basis of the Schw arzschild 
linear programming method ^ISchwarzschildlll979^ . Thus, it 
is not surprising that just for this met hod most signifi- 
cant developments occured in last decade. iRix et all (Il997f) 
and Crotton ct al. ( 1999_) have developed this method in or- 
der to calculate line-of -sight velocity profil e s. Th ereafter, 
ICappellari et alJ ll2002h and IVerolme et"ai] ll2002ll gener- 
alized it for an arbitrary density distribution linking it 
with MGE method. A modification of the least-square algo- 
rithm was done by Krainovic ot al. ( 2005) . Interesting com- 
parisons of the results of the Schwarzschild method with 
phase density calculations within a two-integral app roxi- 
mation have been mad e bv lvan der Marel et alj lll99d) and 
iKrainovic et al.l (|2005"). 

At present, nearly all dynamical models have been ap- 
plied for one-component systems. However, the structure of 
real galaxies is rather complicated - galaxies consist of sev- 
eral stellar populations with difi^erent density distribution 
and different ellipticities. In addition, in different compo- 
nents velocity dispersions or rotation may dominate. 

In our earlier multi-c o mpone nt al models (see 
[ Tenies. Hand fc Einastol Il994 Il998l : lEinasto fc TeniesI 
Il999r) we approximated fiat components with pure rota- 
tion models and spheroidal components with dispersion 
dominating kinematics. For spheroidal components mean 
velocity dispersions were calculated only on the basis of 
virial theorem for multi-componental systems. These mod- 
els fit central velocity dispersions, gas rotation velocities 
and light distribution with self-consistent models. 



In the present paper, we construct a more sophisticated 
self-consistent mass and light distribution model. We de- 
cided to base it on the Jeans equations. For all visible com- 
ponents, both rotation and velocity dispersions are taken 
into account. The velocity dispersion ellipsoid is assumed 
to be triaxial and line-of-sight velocity dispersions are cal- 
culated. Mass distribution of a galaxy is axisymmetric and 
inclination of the galactic plane with respect to the plane of 
the sky is arbitrary. 

In order to discriminate between DM and visible mat- 
ter, it is most complicated to determine the contribution of 
the stellar disc to the galactic mass distribution. Quite of- 
ten the maximum disc approximation is used. In the present 
paper, we attempt to decrease degeneracy, comparing calcu- 
lated models with the observed stellar rotation curve, veloc- 
ity dispersions along the major axis and in addition, along 
several cuts parallel to the major and minor axis. In the case 
of two-integral models for edge-on galaxies this allowed to 
constrain possible dynamical models {Mcrrificld 1991). 

First measurements of veloci ty dispersions along sev- 
eral s lit positions we re made bv [ K ormen dv fc IllingworthI 
Jl982h and llllungworth fc Schechteij Jl982t) . 

Later, simi lar me asurements were perform ed by 

Binnev et al. ('1990"), Fisher. lUingworth fc Franx ('1994, 



Statler fc Smecker-Hane [199 9), Cappellar i et al. (2002). In 
recent years, with the help of integral field spectroscopy, 
complete 3D velocity and dispersion fields have been 
measured already for several tens of galaxies. 

We apply our model for the nearby spiral Sa galaxy 
M 104 (NGC 4594, the Sombrero galaxy). This galaxy is 
suitable for model testing, being a disc galaxy with a sig- 
nificant spheroidal component. The galaxy has a detailed 
surface brightness distribution and a well-determined stel- 
lar rotation curve. M 104 has a significant globular cluster 
(GC) subsystem. And as most important in our case, the 
line-of-sight velocity dispersion has been measured along the 
slit at different positions parallel and perpendicular to the 
projected major axis. 

We construct the model in two stages. First, a surface 
brightness distribution model is calculated. Here we distin- 
guish stellar populations and calculate their structure pa- 
rameters with the exception of masses. In the second stage, 
we calculate line-of-sight velocity dispersions and the stellar 
rotational curve and derive a mass distribution model. 

Sections 2 and 3 describe the observational data used in 
the modelling process and construction of the preliminary 
model. In Section 4 we present the line-of-sight dispersion 
modelling process. Section 5 is devoted to the final M 104 
modelling process. In Section 6 the discussion of the model 
is presented. 

Throughout this paper all luminosities and colour in- 
dices have b een corrected for absorptio n in o ur Galaxy 
according to ISchleeel. Finkbeier fc Davi3 (119981) . The dis- 
tance to M 104 has been taken 9.1 Mpc, correspond- 
ing to the scale 1 arcscc — 0.044 kp c ( Ford et al.l Il99^ : 
iLarsen. Forbes fc Brodi6..20ol iTonrv et al. .200]]) . The an- 
gle of inclination has been taken 84°. 
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Table 1. Photomctrical data 



Reference 


Faintest 


Colour 




isophote 


system 


Soinrad et al. (197H) 




R 


Burkhead dgrQ) 


28.2 


B 


Boroson fl981^ 


24.7 


B 


Hamabe & Okamura fl982~) 


22.7 


B 


Beck et al. f 1984') 


27.3 


B 


^^^'vis Sz FrccrQSyii ^^^^^^ 




V 


Burkhead (1986) 


23.1 


BVRI 


Kent (1988') 


23.3 


R 


Kormendv ('1988') 


19.2 


V 


Crane et al. C 1993') 


17.4 


B 


Emsellem et al. ('1994') 


18.5 


R 


Emsellem et al. (1996) 


17.8 


V 



2 OBSERVATIONAL DATA USED 

By now a surface photometry of M 104 is available in 
UBVRIJHK colours. In the present study, we do not use 
the [/-profile, as this profile has a rather limited spatial ex- 
tent and is probably most significantly distorted by absorp- 
tion. In certain regions also the B-profile is probably in- 
fluenced by absorption, but the B-profile has the largest 
spatial extent and we decided to use it with some cau- 
tion outside prominent dust lane absorption. The JHK- 
profiles have a rather limited spatial extent and resolu- 
tion and we decided not to use them in here. In this way 
the surface brightness profiles in BVRI colours were com- 
piled. Different colour profiles help to distinguish stellar 
populations and allow to calculate corresponding M/L ra- 
tios, and thereafter, colour indices of the components. Ta- 
ble presents references, the faintest observed isophotes 
(mag arcsec"^), and the colour system used. Different R 
colour system profiles are tra nsferred into the Co usins sys- 
tem, using th e calibration b vl^ei fc GunnI 1^9^. The ob- 
servations by ISpinrad et alJ (119781) were made without ab- 
solute calibration. They w ere calibrated with th e help of 
other R colour observations. iHes fc Peletieil (Il993l) observed 
M 104 in BVRI colours but in their paper only colour in- 
dices are given and we cannot use them in here. 

The composite surface brightness profiles in the BVRI 
colours along the major and/or the minor axes were derived 
by averaging the results of different authors. Due to dust 
absorption lane, surface brightnesses only on one side along 
the minor axis have been taken into account. All the surface 
brightness profiles obtained in this way belong to the ini- 
tial data set of our model construction. To spare space, we 
present here the surface brightness distributions in B and 
R only (Fig. Q upper panels), and the axial ratios (the ratio 
of the minor axis to the major axis of an isophote) (Fig. 
lower panel) as functions of the galactocentric distance. 

The observed surf ace density distr i bution of GC candi- 
dates was derived b y Bridges fc Hanej ^ll992^ . ILarsen et alJ 
(|200ll) and lRhode fc Zepi Ii2004f) . The derived distributions 
were averaged, taking into account different background lev- 
els. The resulting surface density distribution of GC candi- 
dates is given in Fig. |5| by filled circles and was used to 
constrain stellar metal-poor halo parameters. Line-of-sight 




0.01 0.1 1 10 100 




26 I ■ i ' ' 

0.01 0.1 1 10 100 




0.0 I 

0.01 0.1 1 10 100 

R (kpc) 

Figure 1. Upper panels: the averaged surface brightness profiles 
of M 104 in the B and R colours. Filled circles - observations, 
solid line - model, dashed lines — models for components (n - the 
nucleus, b - the bulge, h - the metal-poor halo, d - the disc). 
Lower panel: the axial ratios of M 104 isophotes as a function of 
the galactocentric distance. 



velocities of GCs were measured bv lBridges" et al.l II1997D and 
the calculated mean velocity dispersion of GC subsystem 
(tgc = 255 km s~^ was derived. 

Rotation velocities of stars and line-of-sight velocity dis- 
persion profile along the major axis in very good seeing con- 
dit ions (0.2-0.4 arcsec) for the central regions was obtained 
bv iKormendv et all lll99d) with HST an d CFHT. In addi- 
tion, t h e central regions were m easured bv lCarter fc Jenkins! 
|199^ . lEmsellem et alj (ll99(J) . In the central and inter- 
mediate distance interval , dispersions and stellar r otatio n 
have been me a sured by iKormendy fc Illingwor tb| (|l982l) . 
iHes fc PeletTeil (Il993() and Ivan der Marel et al.l ^ll994^ ■ We 
averaged the stellar rotation velocities at various distance 
intervals with weights depending on seeing conditions and 
velocity resolution, and derived the stellar rotation curve 
presented by filled circles in Fig. |H] Averaged in the same 
way line-of-sight velocity dispersions along the major axis 
are presented by filled circl es in Fig. B] 

Tri addition, iKormendv fc Illingworthl lll982h derived 
dispersion profiles along several slit positions (at 0, 30, 40, 
50 arcsec parallel and at 0, 50 arcsec perpendicular to the 
major axis) in the bulge component. 'We use them in mass 
distribution model calculations (filled circles in Figs. llOllllll . 

Ionized gas radial velocities were obtained and the 
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Figure 2. The surface density distribution of GCs in M 104. Tlie 
observations are presented by filled circles. Tlie continuous line 
gives our best-fitting model distribution for the halo. 



Figure 4. Observed line-of-sight velocity dispersions of M 104 
along the major axis. Filled circles — observations, solid line - 
model. 




p^l^pi^^ 4 s b Figure 5. Calculated circular velocity for the best-fitting model 

of M 104 (solid line). Dashed lines give circular velocities for com- 
Figure 3. The rotation curve of M 104 on the basis of stellar ponents (dm - dark matter). Observed gas velocities are given by 

rotation velocities. Filled circles - observations, solid line - model, fiUed circles. Only outermost points are given where stellar mo- 

dashed lines - models for components. tions are not known. 



rotation curve wa s constructed by ISchweizeil il978h and 
iRubin etim lll985l) . HI velo cities at 11.4 arcsec (0.5 kpc) 
resolution were obtained by iBaiaia et all ^si). Unfortu- 
nately, we have no detailed information about the gas ve- 
locity dispersions. When using gas rotation velocities, often 
an assumption is made that gas dispersions are small when 
compared with rotation velocities, and in this way, rotation 
velocities are taken to be circular velocities. However, in the 
case of M 104, up to distances ~ 3 kpc, rotation velocities 
of stars and gas are comparable and thus we may expect 
also dispersions to be comparable, and therefore, gas dis- 
persions cannot be neglected. For this reason, we cannot 
use gas rotation velocities directly in fitting the model. We 
use gas rotation only to have an approximate mass distribu- 
tion estimate at large galactocentric distances where stellar 
rotation and dispersion data do not extend. In these outer 
parts, the velocities from different studies were averaged and 
the resulting gas rotation velocities are given by filled circles 
in Fig. 1^1 



3 SURFACE BRIGHTNESS DISTRIBUTION 
MODEL OF THE M 104 

To construct a model of the M 104 galaxy, we limit the main 
stellar components to the central nucleus, the bulge, the disc 
and the metal-poor halo. To construct a dynamical model 
in the following sections, a DM component - the dark halo - 
must be added to visible components. 

To construct the hght distribution model, the surface 
luminosity distribution of components is usually approxi- 
mated by the Sersic formula (Scrsic 1968). If, in addition 
to the photometrical data, kinematic data are also used, the 
corresponding dynamical model must be consistent with the 
photometry, i.e. the same density distribution law must be 
used for rotation curve modelling (and for the velocity dis- 
persion curve, if possible). For spherical systems, an expres- 
sion for cir cular velocity with an int eger Sersic index can 
be derived IIMazure fc Capelatdl2002^ . For a non-integer in- 
dex and ellipsoidal surface density distribution, a consistent 
solution for rotation curve calculations is not known. 

In the present paper, the density distribution parame- 
ters are determined by the least squares method and may 
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have any value. In addition, our intention is to use the model 
also for velocity dispersion calculations. 

For the reasons given above, we decided to construct 
models starting from a spatial density distribution law for 
individual components, which allows an easier fitting simul- 
taneously for light distribution and kinematics. 

In such models, the visible part of a galaxy is given as 
a superposition of the nucleus, the bulge, the disc and the 
metal-poor halo. The spatial density distribution of each 
visible component is approximated by an inhomogeneous el- 
lipsoid of rotational symmetry with the constant axial ratio 
q and the density distribution law 

Z(a) =Z(0)exph(a/(fcao))'/'^], (1) 

where 1(0) — hL/{ATvqaQ) is the central density and L is 
the component luminosity; a — + z^/q'^, where R and 

z are two cylindrical coordinates, ao is the harmonic mean 
radius which characterizes rather well the real extent of the 
component, independently of the parameter N. Coefficients 
h and k are normalizing parameters, depending on A^, which 
allows the density behavior to vary with a. The definition of 
the normalizing paramete rs h and k and their calculation is 
described in appendix B of lTenies et alJ il994) . Equation Q 
allows a sufficiently precise numerical integration and has a 
minimum number of free parameters. 

The density distributions for the visible components 
were projected along the line-of-sight, and their sum gives 
us the surface brightness distribution of the model 



li{a)a da 



Qi J (a2-A2)l/2 



(2) 



where A is the major semiaxis of the equidensity ellipse of 
the projected light distribution and Qi are their apparent 
axial ratios = cos^ S + q^ sin^ S. The angle between the 
plane of the galaxy and the plane of the sky is denoted by S. 
The summation index i designates four visible components. 

For the nucleus and the stellar metal-poor halo, param- 
eters q, ao, N were determined independently of other sub- 
systems. For the nucleus, these parameters were determined 
on the basis of the central light distribution; for the metal- 
poor halo, these parameters were determined on the basis of 
the GC distribution. In subsequent fitting processes, these 
parameters were kept fixed. This step allows to reduce the 
number of free parameters in the approximation process. 

The model parameters q, ao, L, and A*' for the bulge and 
the disc, and L for the nucleus and the halo were determined 
by a subsequent least-squares approximation process. First, 
we made a crude estimation of the population parameters. 
The purpose of this step is to avoid obviously non-physical 
parameters - relation JJJ is non-linear and fitting of the 
model to observations is not a straightforward procedure. 
Next, a mathematically correct solution was found. Details 
of the least squares approximatio n and the general mod - 
elling procedure were described bv lEinasto fc HaudI ijlQSfltl . 
[Tenies et al. (J994, 199^) . 

Total number of free parameters (degrees of freedom) 
in least-squares approximation was 18, the number of ob- 
servational points was 231. Transition from the bulge to the 
disc and from the disc to the metal-poor halo is rather well- 
determined by comparing the light profiles along the major 
and the minor axis (see Fig. lower panel) . Parameters of 



the nucleus are more uncertain because no sufficiently high- 
resolution central luminosity distribution observations are 
available for us. On the other hand our aim is to study gen- 
eral mass distribution in M 104 where nuclear contribution 
in small. For this reason convolution and deconvolution pro- 
cesses were not used in luminosity distribution model and 
in subsequent mass distribution model. 

The final parameters of the model (the axial ratio q, 
the harmonic mean radius ao, the structural parameters N , 
the dimensionless normalizing constants h and k, BVRI- 
luminosities) are given in Table H The model is repre- 
sented by solid lines in Figs. |5| The mean deviation of 
the model from the observations of surface brightnesses is 
(^°b= - ^■"°'^<=') = 0.16 mag. 



4 CALCULATION OF VELOCITY 
DISPERSIONS 

4.1 Basic formulae 

Knowing spatial luminosity densities of the components U (a) 
and ascribing a mass-to-light ratios to each component fi 
{i indexes the nucleus, the bulge, the disc and the stellar 
metal-poor halo), we have spatial mass density distribution 
of a galaxy 



p{a) 



4 

E 



fih{a) + pDM{a) 



(3) 



(pDM(a) is the DM density). On the basis of spatial mass 
density distributions, derivatives of the gravitational poten- 
tial_J£ and ^ can be calculated (see Biimov & Tromaii3 
I1987D . 

In stationary coUisionless stellar systems with axial 
symmetry the Jeans equations in cylindrical coordinates are 



dR 



+ 



dz 



+ P 



R 



dR ' ' 



d{pvRV0) d{pv0V-^) 2p 



d{pVRV^) d{pvi) pVRV 



OR 



dz 



R oz 



where vr, Vz, vg are velocity components. 
The velocity dispersion tensor afj = i 



(4) 
(5) 
(6) 

viv-j in the 



diagonal form for the axisymmetric case can be described by 
four variables: dispersions along the coordinate axis {or, 
and Gg) and an orientation angle a in .Rz- plane (see Fig.EJ. 
Mixed components of the tensor are 



2 f 2 2\ 



where 



1 ~ 2 2cf- 



2 2 



0, 



(7) 



(8) 



As a result of axial symmetry the second Jeans equation 
vanishes. 

Introducing the dispersion ratios 



kz 



kg = 



-I 



(9) 
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Figure 6. Elliptical coordinates {xi, X2) and their relations with 
cylindrical coordinates {R, z) in galactic meridional plane. 



the remaining Jeans equations can be written in a more 



convenient form for us 



dR 



+ 



R 



+ 



Ok 



2 , 9pa| 
po-fl + — 
oz 



dz 



+ 



or) ■ 



9$ 
dR 



dR 



= -P 



dz ' 



where Ve is rotational velocity, 



(10) 



(11) 



(12) 



For each component the rotation velocity have been taken 
Ve = PVc, where Vc is circular velocity and /3 is a constant 
specific for each subsystems. Taking into account the defini- 
tion of the circular velocity we can substitute in Eq. llOt 



R ^ OR' 



(13) 



The Jeans equations IIUH and IIH include unknown 
values kz, kg and 7. Spatial density, gravitational potential 
and rotational velocity can be determined on the basis of the 
galactic surface brightness distribution (Eq. |5J, the Poisson 
equation and the observed stellar rotation curve. Dispersions 
must be calculate d from the Jeans eq uations. 
Following the notation of lLandau fc Livshit j 



cr|j and 



define confocal elliptic coordinates (a;i, X2) 



R' 



where 



+ 



= 1, 



< 4 



1973), we 

as the roots of 

(14) 



(15) 



The third coordinate Xi = 9. Foci of ellipses and hyperbolae 
are determined by ±zo. The relations between elliptic and 
cylindrical coordinates are as follows; 



x\ = \[^+{^'-Az-'z^f''] 

■ [n-{n^-Az^z^of^^] 



2 

^2 



where 

Q = R^ + zq + z^ . 

Later we also need the relation 



tana 



Rz 



(16) 
(17) 

(18) 
(19) 



In this case, the parameter 7 related to the angle be- 
tween the ellipsoid major axis and the galactic disc is 

Rz 



7 ■ 



i?2 



(20) 



The position of foci zq is at present a free parameter, which 
must be determined within the modelling process. For 20 = 
const the orientation of the velocity ellipsoid would be along 
the elliptic coordinates. Velocity dispersions along elliptical 
coordinates (xi, X2, X3) are denoted as (ai, (T2, 0-3). 

In the case of a triaxial velocity ellipsoid, the phase 
density of a stellar system is a function of three integrals 
of motion. For an axisymmetrical system, in addition to en- 
ergy and angular momentum integrals, a third non-cla ssical 
integral is needed. As it was stressed bv lKuzminl ^ll953^ . this 
third integral should be quadratic with respect to velocities 
(in this case minimum number of constraints result for gravi- 
tational potential). On the basis of this assumption . iKuzminI 
(1953) derived a corresponding form of the third integral. 

Sta r ting f rom the form of Kuzmin's third integral, 
lEinastol il970l) derived that dispersion ratios can be writ- 
ten in the form 



-I 



aizp + a2x\ 
aizl -f a2x\ ' 

aizp -I- a2X% 
aizl + a2Z^ + fe2-R2 ■ 



(21) 



(22) 



where ai, 02 and 62 are unknown parameters. As a sim- 
plifying assumption these three parameters were taken in 



pliiymg assum 
lEinastol \W7(t\ 



(23) 



must satisfy the relation llKuzmin„.1961il 

^-1 + ^- 

K12 K13 

This relation was derived by Kuzmin in the case of disclike 
systems and we must keep in mind that therefore our re- 
sults may not be a good approximation far from the galactic 
plane. 

Using the relations between cylindric and elliptic coor- 
dinates we derive 



kg 



' a + k\2 cos a tan a -I- fc 



12 



a + fci2 sin a 



1 + fci2 tan a 
fci3(l + tan^ a) 



' a + fci2 sin a I + fci2 tan a 



(24) 



(25) 



The quantity zq determines the orientation of the veloc- 
ity ellipsoid. For specific density distribution (gravitational 



© 2006 RAS, MNRAS 000. IHO 



Mass distribution of NGC 4 594 7 




Direction to galaxy centre 
A 




2 2 2 I 2-2 

ai = CT^ cos B + sm B. 



Figure 7. Dispersion projection in a plane parallel to galactic 
disc. 
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Figure 8. Projection of dispersions to the line-of-sight. 



potential) forms within the theory of the third integral of 
motion zo = const. In the case of general density distribu- 
tions zq = f{ R,z). For example, it was demonstrated by 
iKuzminI il962tl that in a galactic disc 



R 



(26) 



where pt is total galactic spatial mass density. From Eq. I|26|l 
we can determine Zq 



zl{R,Q) = -R 



4pt{R 



+ R 



(27) 



The dependence of zo on z is derived to have the best-fitting 
with measured dispersions. 



4.2 Line-of-sight dispersions 

Calculated on the basis of hydrodynamic models, dispersions 
cr|j, and ag cannot be compared directly with measure- 
ments. We project the velocity dispersions in two steps (see 
Figs. IZIIHJ- First, we make a projection in a plane parallel 
to the galactic disc. For this we must project aj^ and ag to 
the disc, going along the line-of-sight and being parallel to 
the galactic disc. Projected dispersions are 



2 2^,2 



1 - 



i?2 



(28) 



Second, we must project dispersions and crj to the line-of- 
sight. Designating O as the angle between the line-of-sight 
and the galactic disc, the line-of-sight dispersion af is 



(29) 



To compare the calculated dispersions with the mea- 
sured data, we must calculate averaged along the line-of- 
sight dispersion. Integrating dispersions along the line-of- 
sight we may write 



L{X,Y) 



l{R,z)af{R,z)dl, 



(30) 



where 1{R, z) denotes galactic spatial luminosity density, 
and L{X, Y) is the surface luminosity density profile (please 
note that integration dl means integration along the line-of- 
sight). Changing variables in the integral to have integration 
along the radius, we obtain 



^LtiX, Y) 



where 



R 



L{X,Y) J cos eVR'^ ~ 



* = 1{R, zi)af{R, zi) + 1{R, Z2)af{R, 22), 



dR, (31) 
(32) 

Equation 13111 gives the line-of-sight dispersion for one galac- 
tic component. As in our model we have several components, 
we must sum over all components considering the luminosity 
distribution profile 

1/2 



Zl,2 



± \/r:^-xA tan e. 
V sm 6 / 



(Tobs(X,F) 



J:[U{X,Y) (aL(X,F))J 

i 



(34) 



where i denotes the subsystem and the summation is taken 
over all subsystems. 



5 VELOCITY DISPERSION PROFILES OF 
NGC 4594 

In the present section we apply the above constructed 
model to a concrete galaxy. We selected the Sa galaxy 
NGC 4594 having enough observational data to construct 
a detailed mass distribution model. The velocity dispersions 
for NGC 4594 have been measured also along several slit 
positions outside of the galactic disc. 

To avoid calculation errors, we first made several tests: 
we calculated dispersions for several simple density distribu- 
tion profiles, varied the viewing angle between the disc and 
the line-of-sight, varied density distribution parameters. All 
test results were in accordance with our physical expecta- 
tions. 

In velocity dispersion calculations all the luminosity dis- 
tribution model parameters derived in Section 3 will be han- 
dled as fixed. The visible part of the galaxy is given as a 
superposition of the nucleus, the bulge, the disc and the 
metal-poor halo. The spatial luminosity and mass density 
distributions of each visible component are consistent, i.e. 
their mass density distribution is given by 

p(a) = p(0)exp[-(a/(fea„))^/'^], (35) 

where p(0) = 1{0)M/L = hM / {A-Kqa^) is the central mass 
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Major Axis (kpc) 

Figure 9. Orientation of tiie calculated velocity dispersion ellip- 
soids in galactic meridional plane. 



density and M is the component mass. For designations see 
Eq. JTJ. 

In the case of mass distribution models, a DM compo- 
nent must be added to visible components. The DM distri- 
bution is represented by a spherical isothermal law 



p(a) 



p(0){[l + (^f]-i-[l + (£)2]-i} a 







a > a 



(36) 



Here a is the outer cutoff radius of the isothermal sphere, 
= kaQ. 

Our model includes an additional unknown value - ve- 
locity ellipsoid orientation. We have to find the best solution 
to zq, when fitting the model to the measured dispersions. 

Figure El gives the shape and orientation of the velocity 
dispersion ellipsoid in the galactic meridional [R, z) plane. 

In Figure IIUI calculated line-of-sight dispersions along 
and parallel to galactic major axis are given. In Figure 1111 
calculated line-of-sight dispersions parallel to the minor axis 
are given. 

It is seen that moving further ofi^ from the galactic disc, 
the results become a little different from the data observed. 
One reason may be that we could not find an appropriate 
solution for zq. It is possible to fit the data far from the 
galactic plane with appropriate selection of zq but in this 
case the fit with dispersions along the major axis is not so 
good. As a result, the figures present the best compromise 
solution we could find. 

In addition we must take into account the observed av- 
erage velocity dispersion of GCs (tgc = 255 km s~^. This 
value is clearly higher than the above mentioned last points 
in stellar dispersion curve ~ 160— 180 km s^^ (Fig. llH . All 
these dispersions correspond to a region where DM takes ef- 
fect. For this reason we think that the mean velocity disper- 
sion of GCs and and stellar velocity dispersions far outside 
of galactic plane can be fitted consistently by introducing a 
flattened DM halo density distribution. This kind of mod- 
els have not yet been constructed by us within the present 
algorithm. 

Figure [T^ gives the calculated velocity dispersion in the 
i?z-plane illustrating the behavior of dispersions. 



6 DISCUSSION 

In the present paper we developed an algorithm, allowing to 
construct a self-consistent mass and Ught distribution model 
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Figure 10. Line-of-sight velocity dispersions of NGC 4594 along 
and parallel to major axis. Solid line - calculated model disper- 
sions, filled circles - observations. 



and to calculate projected line-of-sight velocity dispersions 
outside galactic plane. We assume velocity dispersion ellip- 
soids to be triaxial and thus the phase density is a function 
of three integrals of motion. The galactic plane may have 
an arbitrary angle with respect to the plane of the sky. The 
developed algorithm is applied to construct a mass and Ught 
distribution model of the Sa galaxy M 104. In the first stage 
a luminosity distribution model was constructed on the basis 
of the surface brightness distribution. The inclination angle 
of the galaxy is known and the spatial luminosity distribu- 
tion can be calculated directly with deprojection. Using the 
surface brightness distribution in BVRI colours and along 
the major and minor axis, we assume that our components 
represent real stellar populations and determine their main 
structure parameters. In the second stage, the Jeans equa- 
tions are solved and the line-of-sight velocity dispersions 
and the stellar rotation curve are calculated. Observations 
of velocity dispersions outside the apparent galactic major 
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Table 2. Calculated model parameters. 



Popul. 


M 


ao 


9 


N 


k 


h 




Lb 


Lv 


Lr 


Li 


Nucleus" 


0.001 


0.0015 


0.99 


3.0 


0.00297 


314.3 




3.9E-4 








Bulge 


3.4 


0.28 


0.54 


2.1 


0.03539 


44.51 


0.7 


0.34 


0.48 


0.56 


0.92 


Halo 


7.4 


11.0 


0.75 


4.0 


1.465E-4 


2807. 


0.3 


3.2 




4.4 




Disc 


12.0 


3.4 


0.25 


0.78 


0.7429 


2.607 


0.88 


1.6 


2.3 


2.45 


5.0 


Dark matter 


180 


40.0 


1.0 




0.1512 


14.82 













Masses and luminosities are in units of 10^" Mq and 10^" Lq; components radii are in kpc. 
" A point mass IO^'Mq have been added to the center of the galaxy. 




Figure 11. Line-of-sight velocity dispersions of NGC 4594 par- 
allel to the minor axis. Solid line — calculated model dispersions, 
filled circles - observations. 



axis allow to determine the velocity ellipsoid orientation, 
anisotropy and to constrain DM halo parameters. 

The total luminosity of the galaxy M 104 resulting 
from the best-fitting model is Lb = (5.1 ± 0.6) • 10^°Lq, 
Lr = (7.4 ±0.7) • 10^°Lq. The total mass of the visible mat- 
ter is Mvis = (22.9±3.2)-10^°Mq, giving the mean mass-to- 
light ratio of the visible matter M/Lb = 4.5 ± 1.2 M0L~\ 

M/Lr = 3.1 ± 0.7 MqLq. The surface brightness distri- 
butions in V and I have not sufficient extent to determine 
the luminosities of the stellar halo and we do not either give 
galactic total luminosities in these colours and correspond- 
ing M/L ratios. Calculated from the model, the Lb coin- 
cides well with the total abs olute magni t ude M b ~ —21.3 
(= 5.2 • 10^" L0) obtained bv lFord et all (Il99d) . 

In our model, the mass of the disc is Mdisc = 12 ■ 
10^° Mq. This coincides rather well with the disc mass 
11.4 ■ 10^°Mq calculated with the help of To omre's stabil- 
ity criterion bv lvan den Burg fc Shan el lll986ll and wit h the 
mass 9.6 • 10^° Mr:^ derived bv lEmsellem et all On 
the other hand. lEmsellem et al.j lll994i) derived for the bulge 
mass > 5 • IO^'^Mq, giving Mdisc /Mbuige = 0-2. This issimi- 
lar to the value 0.25 derived bv Ijarvis fc FreemanI il985l) . 
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Figure 12. Projected line-of-sight dispersions (in km s ^) in 
galactic meridional plane. 



but this is much less than Mdisc /Mphcr ~ 1-1 resulting 
fro m our model. An explanati on may be that in the mo dels 
bv lEmsellem et all (jl994) and lJarvis fc FreemanI (|^8^ no 
DM halo was included and hence the extended bulge mass 
is higher. 

In our model, the disc is rather thick (q — 0.25). How- 
ever, disc thickness can be easily reduced to q = 0.15 — 0.2 
when taking the galactic inclination angle instead of 5 = 84° 
to be 83-82°. Other parameters remain nearly unchanged. 
At present this was not done. 

Derived in the present model bulge parameters can be 
used to compare them with the results of chemical evolu- 
tion models. Our model gives M/Lv = 7.1 ± 1.4 MqLq^ 
and (B — V) = 1.06 for the bulge. Compari ng spectral line 
intensities with chemical evolution models, IVazde kis et alJ 
(,1997) obtained for the bulge region t he metallicity Z = 
0.03 a nd the age 11 Gyrs. According to lBruzual fc Charlod 



i2003r . these parameters give M/Lv = 7 — 8 MqLq and 
[B -V) = 1.06 - 1.08 for simple stellar population (SSP) 
models. Bulge parameters from our dynamical model agree 
well with these values and suggest that our model is realistic. 

In our calculations, we corrected luminosities from 
the absorption in the Milky Way only and did not take 
into account th e inner absorption in M 104. According to 
lEmsellem et al.l lll996ll . absorption in the centre may be 
at least Ay ~ 0.13 mag and thus M/Ly ~ 6.3 for the 
bulge. This is slightly to o small when compared with the 
iBruzual fc Charlod ll2003l) SSP models. However, decreasing 
the bulge age to 10.5 Gyrs allows to fit the results. 

Rather sophisticated models of M 104 have been 
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constructed by 


Emsellem et al.l 


fEmsellem & Ferruit 




I2OO0I). Due 



(1199^ 

to our 



199?) and 
diflterent ap- 
proaches, it is difficult to compar e our componen t s and 
their parameters with those of lEmsellem et al.l il994 . 
Il996l) . On the basis of the data used by us, we had no 
reason to add an additional inner disc or a bar to the bulge 
region. However, we did not analyse I and H colours and 
ionized gas kinemat i cs in inner regions as it was done by 
lEmsellem fc Ferruid ll200Cfl . Modelling of gas kinematics in 
central regions is beyond the scope of the present paper as 
gas is not collision- free. 

On the basis of velocity dispersion observations only 
along the major axis it is difficult to decide about the pres- 
ence of the DM even when dispersions extend up to 2-3 Rc 
(Samurovic & Danziacr 2005). In the case of M 104, addi- 
tional dispersion measurements can be used. Velocity disper- 
sions in the case of the slit positioned parallel and perpen- 
dicular to the galactic major a xis, have been measured by 
iKormendv fc IllingworthI lll982h . The calculated mass distri- 
bution model describes rather well the observed stellar ro- 
tation curve and line-of-sight velocity dispersions. Only the 
two last measured points at a cut 50 arcsec perpendicular to 
the major axis deviate rather significantly when compared 
to the model. On the other hand, in addition to stellar veloc- 
ity dispersion measurements, the mean line-of-sight velocity 
dispersio n of the GC subsyst em a = 255 km s~^ was mea- 
sured bv {Bridges et alJ lll997l) . This corresponds to GCs at 
average distances 5-10 kpc from the galactic center and is in 
rather good agreement with the dispersions calculated from 
the model. 

In the best-fitting model the DM halo harmonic mean 
radius ao = 40 kpc and M = 1.8 ■ IO^^Mq giving slightly 
falling rotation curve in outer parts of the galaxy (Fig. |^. 
The central density of the DM halo in our model is p(0) = 
0.033 Mqpc~^, being also slightly less than it was derived 
axies (p(0) = 0.012-0.028 M0pc"^ 
On the other hand, the result fits 
with the limits derived bv l^riell o & Salucci (2001) for local 
galaxies p(0) = 0.015 - 0.050 M0pc~^. 

An essential parameter in mass distribution determina- 
tion is the inclination of the velocity disp ersion ellipsoid with 
respect t o the galactic plane (see e.g. iKxiiiken fc Gilmor j 
Il989l : iM^rrifield 1991) . Velocity dispersion ellipsoid inclina- 
tions calculated in the present paper are moderate, being 
^ 30°. In a sense, our approach to the third integral of stel- 
lar motion is similar to that by Kent fc dc Zocuw ( 199ll) ~ 
the local Stackel fit. In their modelling of the local Milky 
Way structure, they derived that at < 2 < 600 pc and 
6.8 < R < 8.8 kpc, the inclination of the velocity dispersion 
ellipsoid is less than z/R, and they studied the correspond- 
ing correction in detail. In our model, in the same distance 
regions (although the Milky Way and M 104 are not very 
similar objects), inclination correction values are slightly 
smaller. Variations of the corrections with R and z are quali- 
tatively similar. The largest difference is the variation of the 
correction value with z, for which Kent & dc Zccuw ( 1991) 
obtained an increase by 0.1, when moving from z = to 
z = 0.6 kpc, but in our model, corresponding increase was 
only by 0.01. In our model, a significant increase of the el- 
lipsoid inclination angle begins at larger z, which can be ex- 
plained by higher thickness of the disc component of M 104 
(g — 0.25). In a following paper we intend to construct sim- 
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Figure 13. Dispersion ratios in galactic meridional plane. Coor- 
dinate z is in kpc. 



ilar models for other galaxies with velocity dispersion mea- 
surements outside the galactic major axis. This may lead 
to more firm conclusions about the inclination of velocity 
dispersion ellipsoids outside the galactic plane. 

Ratios of the line-of-sight velocity dispersions are given 
in Fig. 1131 It is seen that velocity dispersion ellipsoids are 
quite elongated - anisotropies in the symmetry plane at 
outer parts of the galaxy are < 0.5. 

ModeUing the disc Sb galaxy NGC 288 within 
a constant velocity ellipsoid inclin ation approximation, 
iGerssen. Kuiiken fc Merrifieldl 1^93) estimated that the 
disp ersion ratio Uz / o'r = 0.70. Within epicycle approxima- 
tion l^stfeUeFal] ll2005ll derived for Sb galaxy NGC 3982 
the dispersion anisotropy <Jz/o-r ~ 0.73. These two galax- 
ies are morphologically close to the Sa galaxy mod- 
elled in the present paper, and it is seen that disper- 
sion ratios are more anisotropic in our case. In gen- 
eral, a r adially elongated dispersion ellipsoid is rather 
common llShapiro. Gerssen fc van der Marell2003l) . On the 
other hand, there a re exceptions - galaxy NGC 3949 
iWestfall et al.l (|200^ has the dispersion ratio a^joR = 
1.18). 

By using the quadratic programming method iDeionghd 

Il989l) . the distribution function within the three-integral 
approximation has b een numerically calculated for the S O 
galaxy NGC 3115 bv lEmsellem. Deionghe fc Bacon! Jl999tl . 
Assuming some similarity between SO and Sa galaxies, it is 
interesting to compare the derived velocity dispersion be- 
haviour outside the galactic plane. Although detailed com- 
parison is difficult a similar structure of isocurves is seen. 
But the results disagree in values of ffz/gfl. Our ratios 
are radially elongated, the ratios bv lEmsellem et all (^9^ 
are vertically elongated. Qualitatively this is in agreement 
with the general trend that galaxies of earlier morpholog- 
ical type have larger Gz/o-r ratio (Shauiro ct al. 2003). At 
larger R, the dispersion ratio OzjoR decreases. This is in 
agreement with the decrease of dispersion ratios due to the 
decrease of the role of interactions wit h molecular clouds 
at greater galactocentric distances (see Ijenkins fc BinnevI 
[T99O1') . The use of this explanation in the case of gas-poor SO 
galaxies is not clear. At greater distances z both our and 
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lEmsellem et alJ (^9^ dispersion ellipsoids become more 
spherical. 

By using the Schwarzschild method, dispersion ra- 
tios for E5-6 galaxy NGC 3377 have been calculated by 
ICopin. Cretton fc EmsellemI i2004) . Taking into account re- 
lation between spherical and cylindrical coordinates, the be- 
haviour of the dispersion ratios as a function of 7? and z near 
the galact ic plane is in appro ximate accordance (dispersion 
ratios bv ICopin et al.l (|2to3) are slightly more spherical). 
Rotational properties of elliptical and Sa galaxies are too 
different to compare the ratios ag/oR. Unfortunately, it is 
not possible to compare also orientations of velocity disper- 
sion ellipsoids. 
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